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Abstract 

So far all known singularity-free cosmological models are cylindri- 
cally symmetric. Here we present a new family of spherically sym- 
metric non-singular models filled with imperfect fluid and radial heat 
flow, and satisfying the weak and strong energy conditions. For large 
t anisotropy in pressure and heat flux tend to vanish leading to a per- 
fect fluid. There is a free function of time in the model, which can be 
suitably chosen for non-singular behaviour and there exist multiplicity 
of such choices. 
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1 Introduction 



Although the present day observations indicate that the Universe at large 
scale is homogeneous and isotropic and it is well described by the standard 
Friedman-Robertson- Walker (FRW) model. It is also a well recognised view 
that consideration of inhomogeneity at early times is well in order for generic 
initial conditions and for facilitating formation of large scale structures in 
the Universe. Further it is now known for some time [1-4] that inhomoge- 
neous spacetime admits a family of perfect fluid cosmological models without 
the big- bang singularity. Note that these models are exact solutions of the 
Einstein equation satisfying the causality and energy conditions, yet there 
is no divergence of any physical and geometric parameters throughout the 
spacetime. This happens because spacetime does not admit compact trapped 
surfaces which invalidates application of the singularity theorems [5]. 

From the Raychaudhuri equation [6], it is clear that singularity can be 
avoided only if acceleration or rotation is non-zero. In cosmology rotation is 
generally not favoured and hence acceleration must be present if singularity 
is to be avoided. This means spacetime has to be inhomogeneous. In all the 
known non-singular solutions, shear is also non-zero, indicating anisotropy as 
well. It can be proved for a general orthogonal metric separable in space and 
time variables [4-5] and for a G2- symmetric perfect fluid model that presence 
of shear is essential for presence of acceleration [7] . Though shear contributes 
positively to the effective gravitational charge density in the Raychaudhuri 
equation, its dynamical role as making collapse incoherent can combine well 
with acceleration in avoiding singularity. Thus shear also seems to be playing 
very important role in non-singular models. 

So far all known non-singular solutions are inhomogeneous and isotropic as 
well as have cylindrical symmetry [1-4,8-9]. It has perhaps their cylindrical 
symmetry that comes in the way of their application in practical cosmology. 
Since the Universe is known to be spherical to a good degree, it is pertinent 
to ask whether it is possible to have a spherically symmetric cosmological 
model? In this note we give a singularity-free prescription for imperfect fluid 
with heat flux. This is a prescription rather than a solution for no equation 
has been solved. A general spherically symmetric metric has only four inde- 
pendent stresses, which could be interpreted as density, radial and transverse 
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pressures, and radial heat flux. This will be true for any spherically symmet- 
ric metric. The question is to give a prescription which is free of singularity 
and has proper desired behaviour for physical parameters. This has been 
achieved by letting the Tikekar static model [10], which is a particular case 
of the Tolman solution [11], expand. 

Our model represents a spherically symmetric cosmological universe filled 
with imperfect fluid having unequal radial and transverse pressure and ra- 
dial heat flow. The spacetime satisfies the weak and strong energy conditions. 
It is free of any kind of singularity as all the physical as well as kinematic 
parameters remain finite and regular in the entire range of the variables t 
and r. It has the typical behaviour of a non-singular model; density and 
pressure vanishing for large t and r and being maximum at t — and r = 0, 
expansion parameter and radial heat flux change their sense at t — 0, while 
acceleration tends to zero as r — > 0. There is however an unusual feature 
that heat flows radially inward as the universe expands. The heat flow van- 
ishes at r or t — > as well as r or ±t — ► oo and it falls off as t~ A or r -4 
asymptotically. The pressure anisotropy vanishes for r — > as well as for 
r or ±t — > oo. That is, asymptotically it tends to perfect fluid. The shear 
is also non-zero and hence the model is both inhomogeneous and anisotropic. 



2 The model 

The model is described by the spherically symmetric metric, 

ds 2 = (r 2 + P)dt 2 - ^—tL dr 2 - r 2 (d9 2 + sin 2 9dcp 2 ) (2.1) 

H + P 

where P = P(t). The imperfect fluid is represented by the energy- momentum 
tensor [12], 

T ik = (p + p)UiU k - pg ik + Ap[ciC k + ^{g ik - UiU k )} + 2qc {i u k) (2.2) 

where Ui and q are respectively unit timelike and spacelike vectors, p the en- 
ergy density, p the isotropic fluid pressure, Ap the pressure anisotropy and 
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q heat flux. 



We employ the comoving coordinates to write U; L = y^oo^ and take q = 
y/giiSj. The kinematic parameters; expansion, shear and acceleration for 
the metric (1) read as follows: 



_p r 2 



2(2r 2 + P)(r 2 + P) 3 / 2 ' 
Now applying the Einstein equation, we obtain 



2 = \e\ Ur = 



+ P' 



(2.3) 



8vrp = 
8np r - 



2r 2 + 3P 
(2r 2 + P) 2 

1 

" 2r 2 + P 



(2.4) 
(2.5) 



8iip A 



+ 



2r 2 + P 4(2r 2 + P)(r 2 + P) 2 

-Pr 



2P — 



8-rrq 



(9r 2 + 5P)P 2 
(2r 2 + P)(r 2 + P) 



(2r 2 + P) 3 / 2 (r 2 + P)' 
The pressure anisotropy Ap = p r — p± is given by 



(2.6) 
(2.7) 



8vr Ap 



4(2r 2 + P)(r 2 + P)^ 



2P- 



(9r 2 + 5P)P 2 
(2r 2 + P)(r 2 + P) 



(2i 



Now we just need to make a suitable choice for P(t) such that all physical 
and kinematic parameters remain finite and regular. Any even function of 
t without zero would be an appropriate choice, for instance P = a 2 + f3 2 t 2 
or cosh a kt and so on. Thus there exists a family of non-singular models. 
In general P(t) is free similar to the scale factor in FRW models. In here 
there is a great deal of freedom and one may as well ask that any spherically 
symmetric metric could be viewed as representing an imperfect fluid with 
radial heat flux. This is however true. Even then it is a non-trivial matter to 
have acceptable behaviour for model satisfying energy conditions and being 
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free of singularity. As far as we know this is for the first time such a proposal 
has been made. 

Let us consider the simplest case, P = a 2 + f3 2 t 2 . It is easy to verify that all 
the above physical parameters as well as the kinematic parameters remain fi- 
nite and regular for the entire range of variables; — oo < t < oo, < r < oo. 
This indicates clearly that the model is free of singularity. All the above 
parameters tend to zero for t — > ±00 and/or r — > 00. Asymptotically it 
is low density universe, which when contracts, attains the maximum density 
at t — and r = 0, specified by the parameter a(p max = 3/8na 2 ), which 
can be chosen as small as one pleases to have as large p m ax as desired. For 
large t,/3 will correspond to the Hubble parameter. At t = 0, the expansion 

parameter changes sign (contraction < > expansion) and the acceleration 

tends to zero as r — > 0. The universe starts from low density, contracts to 
high density and again expands to low density state without encountering 
singular behaviour of any kind. Anisotropy in pressure Ap tends to zero for 
r — ► as well as asymptotically (t — ► ±00 or r — > 00). For small r and 
large t, it approximates to the radiation universe, p = 3p, while for large r it 
tends to an isothermal stiff fluid with p = p ~ 1/r 2 [13]. p and p r fall off as 
t~ 2 , whereas Ap falls off as t~ 6 and q as t~ 4 , which indicate that fluid turns 
almost perfect (with isotropic pressure) for large enough t. 

We shall now verify that the model satisfies the weak and strong energy 
conditions; T ik ^ l ^ k > and (T ik — \Tg ik )^ l ^ k > for any non-spacelike 
vector £\ We have 

T lk Ci k = n 2 (p + p r - Ap) - m 2 (p r - Ap) 

and 

(2n 2 - m 2 )(p + p r - Ap) + m 2 (2p r - Ap) 

where ^u 1 = n,^ 1 = m 2 and n 2 > m 2 always. It is clear from eqns (4), 
(5) and (8) that the weak and strong energy conditions will be satisfied for 
P(t) = a 2 + (3 2 t 2 . A lengthy but straight forward calculation shows that 
the model cannot however satisfy the dominant energy condition (the vector 



(Tik - 7;Tgik)€£, k — o 
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Tik(, k is non-spacelike). 



All this is very fine but there is a discomforting feature indicated by qd < 
(from (3) and (7)) in general independent of specific choice for P(t). This 
implies that there is a radially inward heat flux for expanding phase and out- 
ward for contracting phase. It may be noted that for the particular model 
under consideration, q = for r or t = and it falls off as r -4 or t -4 for large 
r or t. However, we must confess that this is rather an unusual feature. 

The overall evolution of the model is typical of non- singular models [3,4]; 
asymptotically low density passing through the dense state at t = 0, where 
interchange occurs between expansion and contraction on one hand and be- 
tween inflow and outflow for the radial heat flux on the other. The remarkable 
features of this model are: (a) it is free of big-bang as well as any other singu- 
larity, (b) it is spherically symmetric, which augurs well with the symmetry 
of the realistic Universe, (c) for small r and large t it approximates to the 
radiation Universe, p = 3p, (d) asymptotically (t — ► ±oo and/or r — ► oo), 
the pressure anisotropy and heat flux vanish leading to perfect fluid with 
p = 3p ~ 1/t 2 for large t and p = p ~ 1/r 2 for large r, (e) the parameter a 
defines the maximum density (p m ax = 3/8na 2 for t — 0, r = 0) and (3 will 
correspond to the Hubble parameter for large t, (f) the expansion parameter 
9 and heat flux q change sign at t — and fall off to zero as t — > ±oo and/or 
r — > oo and q6 < 0, (g) the acceleration u r tends to zero as r — > and 
also falls off to zero asymptotically (t — >■ ±oo and/or r — > oo), and (h) 
it is evidently causally stable and obeys the weak and strong energy condi- 
tions. It however does not satisfy the dominant energy condition, which may 
perhaps be responsible for unusual behaviour for heat flow (radially inward 
as the universe expands). It has also been argued that at high densities, the 
possibility of violation of the dominant energy condition cannot be ruled out 
[14]. That means at high density dp /dp may not always measure the true 
signal propogation velocity and hence even if dp/ dp > 1 there is no conflict 
with special relativity. Thus possibility of p > p cannot be ruled out. In the 
context of our model, it suggests that at high density some unusual behaviour 
may be permissible. However the model is free of singularity and satisfies 
the weak and strong energy conditions. 

One of the main objections against the so far known non-singular cosmo- 
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logical solutions was that they did not accord to spherical symmetry and 
hence their application to realistic cosmology was greatly marred. By find- 
ing a family of spherically symmetric non- singular models we have overcome 
this objection quite successfully and hence paving way for their practical ap- 
plication in cosmology. It should be noted that the model is as general as 
FRW with P(t) being free. Non-singular character will however constrain 
the choice for P(t), but even then there would remain good deal of freedom. 
It would therefore be pertinent to examine its cosmological viability. This is 
what we would like to do next. Finally it is quite probable that there could 
similarly exist other families of non-singular spherical models which would 
make it all very exciting. 
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